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We study the dynamics of skyrmions in Dzyaloshinskii-Moriya materials with easy-axis anisotropy. 
An important link between topology and dynamics is established through the construction of un¬ 
ambiguous conservation laws obtained earlier in connection with magnetic bubbles and vortices. 
In particular, we study the motion of a topological skyrmion with skyrmion number Q = 1 and 
a non-topological skyrmionium with Q = 0 under the influence of an applied field gradient. The 
Q = 1 skyrmion undergoes Hall motion perpendicular to the direction of the field gradient with a 
drift velocity proportional to the gradient. In contrast, the non-topological Q = 0 skyrmionium is 
accelerated in the direction of the field gradient, thus exhibiting ordinary Newtonian motion. When 
the applied held is switched off the Q = 1 skyrmion is spontaneously pinned around a fixed guiding 
center, whereas the Q = 0 skyrmionium moves with constant velocity v. We give a systematic 
calculation of a skyrmionium traveling with any constant velocity v that is smaller than a critical 
velocity v c . 


I. INTRODUCTION 

Magnetization structures in the form of topological 
solitons have come to the center of research interest fol¬ 
lowing the fabrication of various ferromagnetic materi¬ 
als. Magnetic bubbles were intensively studied in the 60s 
and 70s for technological applications [T], and theoreti¬ 
cally predicted magnetic vortices were experimentally ob¬ 
served as ground states in mesoscopic magnetic elements 
in the 90s 12113 - Stable topological solitons had been pre¬ 
dicted also in the presence of the Dzyaloshinskii-Moriya 
(DM) interaction jU [5]. They were observed in recent 
years as isolated structures El [7] or forming lattices El- 

IIQ] 

The dynamics of topological magnetic solitons, such as 
bubbles and vortices, has been long recognized to display 
peculiar features El- They are deflected in a direction 
almost perpendicular to an applied magnetic field gradi¬ 
ent pT] and their dynamics is seen to be similar to the 
Hall motion of a charged particle in the presence of a 
magnetic field. A theoretical description links the dy¬ 
namics of topological solitons with their nonzero topo¬ 
logical number, by example of magnetic bubbles mm- 

The opportunity now arises to study peculiar soliton 
dynamics for the case of skyrmions in DM materials. 
These materials are well suited for such studies for two 
main reasons. The DM interaction introduces an intrinsic 
length scale which defines the size of the skyrmion, there¬ 
fore the skyrmion is expected to be robust and remain 
rigid under external probes. In the presence of easy-axis 
anisotropy we have a topological skyrmion with skyrmion 
number Q = 1 as well as non-topological Q = 0 solitons 
(27 r vortices) HU US] thus allowing to explore theoretical 
predictions for dramatically different dynamical behav¬ 
iors. 

Skyrmions could be the stable and robust entities that 
are needed for the technology of recording and transfer¬ 
ring information, currently mainly obtained in magnetic 
media using domain walls. Current induced motion is 


under intensive study and it is a promising technique for 
the manipulation of magnetic information mm- The 
different dynamical behaviors which will be described in 
this paper can help skyrmions emerge as particularly at¬ 
tractive entities for applications in the area of transfer of 
magnetic information. 

Sec. [TT] gives a description of the model, Sec. |III| dis¬ 
cusses static Q = 1 and Q = 0 solutions, Sec. |IV| con- 
tains a review of the theory and gives the main results 
of the dynamics of skyrmions and Sec. [V] contains our 
concluding remarks. 


II. THE MODEL 

We assume a thin film with easy-axis anisotropy per¬ 
pendicular to the rry-plane of the film and with a DM 
energy term (5]. If M = M(x,y) is the magnetization 
vector the energy functional reads 

w = A / d,M ■ dxdy + A J (M 2 + Ml) dxdy 

+ 3 ^ J [( M 3 d 2 - M 2 d 1 )M 3 - (d 2 M 1 - d 1 M 2 )M 3 ] dxdy, 

(1) 

where M s is the saturation magnetization, A is the ex¬ 
change constant, K the anisotropy constant, and D the 
DM constant. Spatial derivatives in Eq. 0 are denoted 
by with n = 1,2 and d\ = d x , d 2 = d y . We have not 
included the energy of the demagnetizing field in Eq. (|T|) 
because it does not affect skyrmion configurations in a 
qualitatively significant way ;22]; it introduces a depen¬ 
dence of the skyrmion size on the film thickness [23] . 
Note that W in Eq. ([I]) is actually the energy per unit 
length along the easy axis perpendicular to the film. 

We define the normalised magnetization m = M/M s , 
so that m 2 = 1. We further use £d = 2A/|Z)| as the unit 
of length, hence the energy measured in units of 2 A is 


2 


I dxdy 


given by (see also [21] ) 

W = ^ J d^m ■ d^m dxdy + ^ J (m 2 + 

+ A J [(toi9 2 - rn 2 di)m 3 - {d 2 m 1 - dim 2 )m 3 ] dxdy, 

( 2 ) 


where 

a = M = ±1 

will be referred to as the chirality and 


K 


D 2 


w, K n = w 


K n 


4 A 


(3) 


(4) 


is the rationalized (dimensionless) anisotropy constant. 
Unless otherwise stated, we choose chirality A = 1 in 
all of our numerical calculations, while k is taken to be 
positive (easy-axis anisotropy). 

The conservative Landau-Lifshitz (LL) equation asso¬ 
ciated with the energy ([2]) is 

dm t 

-at=- mxf ’ (5) 

- SW 

/ = = Am 

dm 


The time variable t is measured in units of tq = 
2AM s /( r yD 2 ) where 7 is the gyromagnetic ratio. 

We recall at this point that magnetic configurations 
are characterised by the skyrmion number defined as 

Q= 4^ 


qdxdy, q = \e^ v m ■ {d v m x d^m), (6) 


where q is called the topological density. The skyrmion 
number Q is integer-valued (Q = 0,±1,±2,...) for all 
magnetic configurations such that m = (0,0, ±1) at spa¬ 
tial infinity. For definiteness we assume m = (0, 0,1) in 
all our calculations. 


We also construct a tensor tr Ml/ defined from 




. ~ SW 

= -/ • d M m = — • d ;l m. 


A formal calculation gives the tensor components 

1 h\j 

(Jn =- (dim • dim — d 2 m ■ d 2 m ) + — (m 2 + m\) 

+ X(mid 2 m. 3 - m 3 d 2 mi) 

(j 12 = — dim • d 2 m + \{m 3 dimi — midim 3 ) 
cr 2 i = - dim • d 2 m + A(m 2 d 2 m 3 — m 3 d 2 m 2 ) 


( 8 ) 


1 fa 

c 22 =- (d 2 m • d 2 m — dim • dim) + — (m 
+ X(m 3 dim 2 — m 2 dim 3 ). 


+ m 2 2 ) 


The topological density q together with the tensor cr M!/ 
provide important theoretical tools for the analysis of 
both static and dynamical properties of the Landau- 
Lifshitz equation mug. 


III. STATIC SKYRMIONS 

The uniform (ferromagnetic) states m = (0, 0 ± 1) are 
the simplest static solutions of the LL equation (|5j) with 
total energy W = 0. For large rationalised anisotropy k 
the uniform state is the ground state of the system. How¬ 
ever, for sufficiently low anisotropy a spiral state with en¬ 
ergy W < 0 becomes the ground state [5] . The transition 
happens at the critical value 


k c = — « 2.4674. 
4 


(9) 


The period of the spiral increases for increasing 
anisotropy and goes to infinity at the critical value k c . 

In order to find nontrivial static solutions which satisfy 
the equation m x / = 0 we may employ a relaxation 
algorithm based on 


dm 

— = - m x (m x /), 


( 10 ) 


nm 3 e 3 

- 2A [d 2 m 3 e± - di m 3 e 2 + (di m 2 - d 2 mi) e 3 ]. 


typically implemented on a 300 x 300 lattice with uni¬ 
form spacing Ax = Ay = 0.1 and Neumann boundary 
conditions. The numerical lattice is large enough so that 
it can be assumed that an infinite film is simulated in 
all presented calculations. An initial spin configuration 


will evolve under Eq. ( 
decreases monotonical 


10| ) in such a way that its energy 
.y and eventually converges to a 
static solution of Eq. ([5]) . 

We now focus on axially symmetric skyrmion configu¬ 
rations. These are conveniently described in terms of the 
standard spherical parametrisation for the magnetization 
given by 

mi = sin 0 cos <&, m 2 = sin 0 sin $, m 3 = cos 0 

(n) 

using the ansatz 


0 = 9(p), $ = </> +7t/2, 


( 12 ) 


where (p, <f>) are polar coordinates. Stationary solutions 
of the energy functional (IT]) then satisfy the ordinary dif- 


(7) 

ferential equation HI El 


d 2 0 1 de f 1 A „ . „ 2A . o /> 

—z -\ ---( k H—- I cos 6 sm# -f-sm 2 9 = 0 (13) 

dp z p dp \ p z J p 

while the skyrmion number defined from Eq. |6]) reduces 
to 

1 POO 1 -1 

q =2 J q -^dp=-[m 3 (oo)-m 3 m ( 14 ) 

where m 3 = cos 9 = m 3 (p) is the third component of 
magnetization. Thus, if Eq. (141 is solved with boundary 
conditions 6(p = 0) = tt and 0[p —► 00) = 0, it leads to a 
static skyrmion with Q = 1. 

Boundary-value problems of the above nature are typ¬ 
ically solved by some sort of a shooting method H2- 
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FIG. 1: The static axially symmetric (Q = 1) skyrmion repre¬ 
sented through the projection (mi, m 2 ) of the magnetization 
vector on the plane for anisotropy k = 3. 


FIG. 3: The static axially symmetric skyrmionium (Q = 0) 
represented through the projection (mi, m 2 ) of the magneti¬ 
zation vector on the plane for anisotropy k = 3. 



P 

FIG. 2: The profiles m 3 = cos# = ms(p) of skyrmions 
(Q = 1) for three values of the anisotropy parameter k = 
2.6, 2.8, 3.0. 


Here we employ a more general (actually simpler) method 
based on the fully-dissipative algorithm (101 which does 
not a priori assume axial symmetry. Instead, we may ini¬ 
tialize Eq. (101 with an essentially arbitrary spin configu¬ 
ration with Q = 1 which then converges to a static Q = 1 
skyrmion that is a local minimum of the energy func¬ 
tional ([ 2 ]). The result is shown in Fig. [I] for anisotropy 
k = 3 and its axial symmetry is evident. Also note that 
the chirality of the skyrmion (counter-clockwise rotation 
in Fig. |T|) is fixed by the sign of the DM constant (A = 1). 
If we choose A = — 1 then the profile m 3 = ms(p) of the 
skyrmion remains unchanged but its chirality is reversed. 
To be sure, the skyrmion number is Q = 1 for either 
choice of chirality. 

Fig- 1 shows the skyrmion profiles for three values of 
the anisotropy constant k. As k approaches k c from 


above, the radius of the skyrmion diverges to infinity. 
There is no evidence for the existence of skyrmions when 
k < K c . 

Other axially symmetric static skyrmions can be found 
if we assume multiple rotations of the magnetization as 
we move radially from the skyrmion center. These were 
first identified in Ref. [14] solving Eq. (13) with bound¬ 
ary conditions 0(p = 0) = kir,9(p —> 00) = 0, with 
k = 1, 2,3,..., where k = 1 is the fundamental (Q = 1) 
skyrmion already discussed, k = 2 leads to a Q = 0 con¬ 
figuration, in view of Eq. (14), and so on. 


In order to apply our relaxation algorithm for k = 2 , 
we first construct a suitable ansatz to be used as initial 
condition in Eq. ( |T()| . We recall the axially symmetric 
<2 = 1 configuration constructed above, which we now 
denote by n = (771,77.2,773). We then apply the transfor¬ 
mation [ 25 ] 

mi = 277.3711, m2 = 2113122, m3 = 211,3 — 1 . ( 15 ) 

The resulting configuration m remains axially symmet¬ 
ric but its skyrmion number is Q = 0 . This follows from 
the fact that m3(p = 0 ) = 1 = 7773 (p —7 00). The mag¬ 
netization rotates to m3 = — 1 at some intermediate ra¬ 


dius. We apply algorithm (|T0) inserting configuration 

gorithm converges to an 


(|15j) as initial condition. The al; 
axially symmetric configuration with <2 = 0, shown in 
Fig. m for k = 3. Such a configuration may be called 
a “skyrmionium” [2'6j because it consists internally of a 
skyrmion and an antiskyrmion. A comparison between 
the skyrmion and the skyrmionium is given through their 
profiles 7773 (p) plotted in Fig. [ 4 ] for the same anisotropy 
parameter k = 3. 

The existence of a static chiral skyrmion has been rig¬ 
orously established within a related model only for <2 = 1 
]27l . We have searched numerically for static solutions of 
E q : § with Q = -1 or 2, which lack axial symmetry, 
but our relaxation algorithm (101 did not converge to 
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FIG. 4: The profiles m 3 = cos# = m 3 (p) of a skyrmion (Q = 
1) and a skyrmionium (Q = 0) for anisotropy k = 3. 


any such static solutions. We will continue this study 
focusing on the axially symmetric configurations shown 
in Figs. Q [3] 

For static solutions such as those discussed in the 
present section the tensor of Eq. © satisfies d\a^\ = 0. 
This yields x^dxcr^x = 0, while integration of both sides 
over the entire plane and an elementary application of 
the divergence theorem leads to [21] 


0) skyrmionium. We shall find that significant differences 
arise in the two cases, as anticipated by an important link 
between topology and dynamics established in our early 
work mm] through the construction of unambiguous 
conservation laws. 

The main result is established rigorously by examin¬ 
ing the time evolution of the topological density q of 
Eq-©- A straightforward application of the LL equation 
© leads to 

Q = • d v m) = e ^ d fl d x a u x (19) 


where the overdot denotes time derivative and ova is the 
tensor defined in Eq. ©. An immediate consequence of 
Eq. (19) is that the integrated topological density, the 
skyrmion number Q of Eq. © is conserved, as expected. 
Furthermore, the appearance of the double spatial deriva¬ 


tive in the right-hand side of Eq. (19) suggests that some 


of the low moments of the topological density are also 
conserved. The lowest nontrivial moments are given by 


1 /j = J x^qdxdy, /x = 1, 2 


( 20 ) 


and their sconservation ( I u = 0) is demonstrated by a 
simple application of Eq. ( [To] ) and the divergence theo¬ 
rem. In order to reveal the physical content of moments 
( f2(~)] ) we note that under a rigid translation of spatial coor¬ 
dinates by a constant vector ( x^ -A a^+c^) the moments 
transform as 


<J[iv dxdy = 0 


(16) 


where the indices /i and v take the values 1 or 2 in any 
combination and thus lead to four independent virial re¬ 
lations that must be satisfied by any static solution. A 
special combination of these relations yields the result 


J (oil + 0 - 22 ) dxdy = 0 => 2W a + 1T D m = 0, 


(17) 


which was obtained through a scaling argument in 
Ref. [29] in analogy to Derrick-Hobart relations [30] |3T|. 
Here, W a is the positive anisotropy energy while Wdm is 
the Dzyaloshinskii-Moriya energy that may be positive 
or negative and thus does not exclude nontrivial static 
solutions. Another special case of Eqs. (16) gives 


Ifi + 47tQc 


'£*• 


( 21 ) 


which indicate an important difference between topolog¬ 
ical (Q ^ 0) and non-topological (Q = 0) magnetic soli- 
tons. The two cases are studied in turn in subsections 
ITVAl and lTVRl 

In order to probe skyrmion dynamics we consider the 
effect of an applied magnetic field h = h(x,y,t) which 
may be a nontrivial function of both spatial and time 
variables. The question is then to predict the behavior 
of the magnetic configuration after the field is turned on. 
The Landau-Lifshitz equation © is then modified by the 
simple replacement f —> f + h. In particular, relation 
(19) becomes 


q = d^dxvvx ~ ■ d v m). 


( 22 ) 


J(cr 12 - a 2 i)dxdy = 0 


J [m 3 (V • m) — (m ■ V)m 3 ] dxdy = 0, (18) 

where V = ei<9i + e 2 d 2 , and it is entirely due to the 
anisotropic DM interaction. 


IV. SKYRMION DYNAMICS 

We now study the dynamics associated with the two 
static solutions calculated in Sec. |III| namely, the topo¬ 
logical (Q = 1) skyrmion and the non-topological (Q = 


As a result the skyrmion number Q is still conserved, as 
expected, but the moments are no longer conserved 
and satisfy 

h- d v m dxdy. (23) 


Yet studying the degree to which the conservation is vi¬ 
olated will give important information on the motion of 
skyrmions in the presence of the applied field. 

The rationalized field h is related to the physical field 
H by 


h=* 

H 0 


D 2 

2 AhqM s 


( 24 ) 
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and we may also write Hq = M s (£ 6 x/^d) 2 where £ ex = 
y/2A/(n>M*) is the exchange length, the most com¬ 
monly used length scale in micromagnetics. 


A. Hall motion of Q = 1 skyrmion 


In view of Eq. (21) the physical interpretation of mo¬ 
ments in Eq. (20) depends crucially on the skyrmion 


number Q. For Q ^ 0, the normalized moments 


R = ^ 1 

M 47tQ 47rQ 


x^qdxdy, 


(25) 


provide a measure of position. The 2D vector R = 
( R x ,R y ) will be referred to as the guiding center of the 
magnetic configuration in question and is conserved in 
the absence of a magnetic field gradient. Thus, a mag¬ 
netic soliton with Q / 0 cannot move freely and is 
spontaneously pinned within the ferromagnetic medium. 
However, motion is possible in the presence of an applied 
magnetic field gradient. 

We consider an initially static skyrmion such as the 
Q = 1 skyrmion shown in Fig. [l] and apply an applied 
magnetic field with a gradient in the x direction: 


ft. = ( 0 , 0, ft), h = gxe x (26) 


where g is the strength of the gradient and a is a constant. 
The field gradient is almost uniform for x < a and it fades 
out for x > a. A pure gradient (ft = gx) may be achieved 
in the formal limit a —> oo. However, such a limit should 
be taken with caution because the field would then reach 
high values at large x and might destroy the uniform 
ground state m = ( 0 , 0 , 1 ) in the presence of dissipation. 

The velocity of the skyrmion guiding center is obtained 
by inserting Eq. (26) in Eq. (23) and applying a partial 
integration 


R-r = 0, 


R =-J_ 

y 47 tQ 


d x h (1 — m 3 ) dxdy. (27) 


The guiding center drifts along the y axis, thus in a di¬ 
rection perpendicular to the applied field gradient. For 
a pure gradient ft = gx , applying Eq. (27) we find the 
simple formula 


Rx = 0 , 


R - 9 V 
v ~ 4 ttQ’ 


(28) 


where g = f (1 — m 3 ) dxdy is the total magnetization. 
For the static Q = 1 skyrmion of Fig. [l] the numerical 
calculation gives g = 15.1, thus the expected velocity is 
R y = —I. 2 O 5 . The above prediction is rigorously correct 
during the initial stages of the process, but deviations 
from a constant velocity are possible at later stages be¬ 
cause the moment g is not by itself conserved. 

The preceding results are verified by a direct numerical 
simulation. We apply a field gradient (26) with g = —0.1 


and a = 10 and solve the initial value problem with the 



FIG. 5: Contour plots for m 3 for a skyrmion under the field- 
gradient (26 1 with g = —0.1, a = 10. (Left) The static Q = 1 
skyrmion of Fig. |T] is placed at the origin at initial time t = 0 . 
(Right) The skyrmion at simulation time t = 30. It moves 
along the y axis perpendicular to the field gradient. The 
contour levels plotted are m 3 = 0.9,0.6,0.3,0.0 (solid lines) 
and m 3 = —0.3,—0.6,—0.9 (dashed lines). 



FIG. 6 : (Left) The trajectory of a Q = 1 skyrmion under the 
influence of the field-gradient (261 with g = — 0 . 1 , a = 10 . 
The guiding center propagates along the y-axis; thus in a di¬ 
rection perpendicular to the applied field-gradient (straight 
black line). The trajectory for the (A', Y) of Eq. ( |29| is the 
cycloid shown by a solid red line. After the applied field is 
suddenly switched off the guiding center ceases to move fur¬ 
ther and the skyrmion trajectory organizes itself in a rota¬ 
tional cyclotron-type motion around the fixed guiding center. 
(Right) A magnification of the latest stages of the process is 
shown. 


initial condition provided by a static Q = 1 skyrmion 
(with k = 3) placed at the origin (0,0). The skyrmion 
moves coherently in the y direction, that is, perpendic¬ 
ular to the field gradient. Fig. [5] shows the skyrmion at 
time t = 0 and at t = 30 when it has drifted approx¬ 
imately Ay = 4 units along the y axis. The guiding 
center, shown in Fig. [ 6 ] (solid black line), moves along 
the y axis with an almost constant velocity. Its initial 
(t < 4) velocity is R y = 0.12 in excellent agreement with 
the prediction obtained from Eq. (28). At later times 
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the guiding center velocity increases to R y = 0.13 while 
a corresponding increase in the moment y is observed. 

One could heuristically invent more quantities to de¬ 
scribe the skyrmion position, whose usefulness depends 
on their connection to actual measurements. A plausible 
measure of position is given by the moments of the third 
component of the magnetization: 

= f x(l -m 3 )dxdy = / y( 1 - m 3 ) dxdy 
f (1 — m 3 ) dxdy ’ j{\ — mz)dxdy 



Fig. § also shows the trajectory (A, Y) obtained through 
the numerical simulation. The position (A, Y ) closely fol¬ 
lows the guiding center, but is also decorated by Larmor 
oscillations. The picture is strongly reminiscent of the cy¬ 
cloidal Hall motion of an electric charge that moves under 
the influence of an electric field and a magnetic field per¬ 
pendicular to the plane of the motion (see Ref. [32] pp 

57 )- 

The preceding analogy with Hall motion is further sub¬ 
stantiated by suddenly switching off the applied mag¬ 
netic field gradient after the skyrmion has moved a dis¬ 
tance, say, Ay = 4 along the y-axis. The guiding center 
then ceases to move further, while the skyrmion orga¬ 
nizes itself in a cyclotron-type rotational motion around 
the pinned guiding center (see Fig. | 6 j) . 


B. Newtonian motion of Q = 0 skyrmionium 


According to Eq. (21 1 the moments are invariant 


under rigid translations for a Q = 0 skyrmion. Therefore, 
unlike the case of the previous subsection for Q ^ 0, the 
conservation laws for I „ do not exclude the possibility 
of free translational motion when Q = 0. Furthermore, 
the linear momentum is defined via the conserved as 

[huh] 


FIG. 7: Contour plot for m 3 for the static skyrmionium (left) 
and the propagating skyrmionium with velocity v = 0.07 
(right). The contour levels plotted with solid lines are m 3 = 
0.9, 0.6, 0.3, 0.0 and with dashed lines m 3 = —0.3, —0.6, —0.9. 


relations generalizing the virial theorems (161 satisfied by 
static solutions. 


An immediate consequence of Eq. (32) is obtained by 


integrating both sides over the entire zy-plane: 

e v \V\Q = 0 (33) 

where Q is the skyrmion number of the magnetic soli- 
ton in question. Therefore, we conclude that topological 
solitons (Q 7 ^ 0 ) cannot be found in rigid translational 
motion ( v\ = 0). This elementary result is in agreement 
with a similar conclusion reached for (Q = 1) skyrmions 
on the basis of the conservation of the guiding center R 
in the absence of an applied field gradient. 

On the other hand, Eq. (33) does not exclude a rigidly 


moving non-topological soliton (Q = 0 ) because it is 
then trivially satisfied for any velocity v\. In this case, 
Eq. ([32]) can be further iterated by multiplying both sides 


with x y and then integrating over the entire xy-plane to 
obtain the virial relations [28] 


Pv = tuvlv, /q v = 1 or 2. (30) 


( v\Px)5 flv - P^Vv = 


cr^ dxdy 


(34) 


The invariance of linear momentum under rigid transla¬ 
tions is in accordance to the properties of ordinary mo¬ 
mentum of a point particle in Newtonian dynamics. 

Newtonian dynamics allows for a steady-state with 
constant velocity when no forces are applied. We thus 
look for steady-states propagating rigidly with constant 
velocity v = (vi,v 2 ) he., m = m{x — vt;v) satisfying 
m = —v\d\m and thus 

v\d\m = mx f. (31) 

Now, take the cross product of both sides with m, then 
the dot product with d„m, to obtain 

e-u\V\q = d\a v \ (32) 


where P y is the linear momentum defined in Eq. (30) 


while indices /i and v take the values 1 or 2 in any combi¬ 
nation and thus lead to four independent virial relations. 
These must be satisfied by any non-topological (Q = 0) 
magnetic soliton moving rigidly with constant velocity 
V = (v 1 ,v 2 ). 


Actual solutions of Eq. (311 are obtained numerically 


by a generalization of the relaxation algorithm ( 10 ) used 
for the calculation of static solutions in Sec. Tml For 
definiteness we assume rigid motion along the cc-axis 
(ui = v,v 2 = 0). Eq. (10) is then generalized accord¬ 
ing to 


dm 

1h 


= —m x ( m x / — v 


dm 

dx 


v = u— P (35) 


where q is the topological density defined in Eq. © and 
a u \ the tensor defined from Eq. Q. Eq. ( f32j ) provides 
the basis for the derivation of a series of interesting virial 


where the velocity v is determined self-consistently in 
terms of an arbitrary input parameter u and the mo¬ 
mentum integral P = P x = I 2 calculated from Eq. (30) 
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FIG. 8: Contour plot for the topological density q of the static 
skyrmionium (left) and the propagating skyrmionium with ve¬ 
locity v = 0.07 (right). The contour levels plotted are chosen 
arbitrarily. Solid lines mean q > 0 and dashed lines q < 0. 


for each time step. It should be noted that the momen¬ 
tum is not conserved by the dissipative dynamics (351, 
hence the “velocity” v evolves together with the spin con¬ 
figuration until they both reach definite terminal values 
that are a local minimum of the Lyapunov functional 
F = W + I(w — P ) 2 . In view of the constraint m 2 = 1, 
which is compatible with Eq. (351, the terminal state will 
satisfy the differential equation 


dm 

dx 


mx f 


(36) 


and thus describes a magnetic soliton that moves rigidly 
along the x-axis with constant velocity v which depends 
on the input parameter u and may be varied accordingly. 

The algorithm is initiated with the static skyrmionium 
of Fig. [3] and a nonzero value for the input parameter u. 
The algorithm converges to a steady-state configuration 
with velocity in the range 


0 < v < v c , v c ~ 0 . 102 . 


(37) 


Fig. 0 shows the static together with a propagating 
skyrmionium via contour plots for the magnetization 
component m 3 . Note that a skyrmionium with v 7 ^ 0 
is no longer axially symmetric. 

A further representation of the skyrmionium is given 
via the topological density in Fig. [ 8 ] The negative lump 
of topological density q moves off the center of the config¬ 
uration for v > 0. This gives manifestly a nonzero value 
P > 0 for the x component of the linear momentum 
(301. For larger velocities the large axis of the elliptically 


shaped contours increases and apparently diverges to in¬ 
finity for v —> v c . We were able to numerically calculate 
the propagating configurations up to v ~ 0 . 102 . 

The energy IT of a propagating skyrmionium can be 
calculated from Eq. (2j) and its linear momentum P = 
( P , 0) is given by Eq. (301. Fig.[9]shows that both W and 
P increase with velocity v. For the static skyrmionium 
the energy W has a nonzero value while P = 0. Both 
quantities diverge to infinity as v —> v c . We could verify 


FIG. 9: (Left) Solid lines show the energy W of Eq. § and 
linear momentum P of Eq. (301 for a steady-state propagating 
skyrmionium as a function of its velocity v. The open circles 
show results from a direct simulation discussed in the text. 
(Right) The energy W versus momentum P for a steady- 
state propagating skyrmionium. The curve is parabolic for 
low momenta in accordance with Eq. (411 but linear for large 
momenta in agreement with Eq. (42 I. 


numerically that the group velocity relation 


dW 

Ip 


= V 


(38) 


holds to a very good accuracy (~ 1 %). 

It is particularly interesting that P is linear for values 
of the velocity v < 0.06, and this motivates the defini¬ 
tion of a mass m for the skyrmionium from the Galilean 
relation 


P = mv 


(39) 


with the value m = 117 extracted from the numerical 
data. Interestingly, the energy at low velocity is given by 


1 

W = Wn + -mv 


(40) 


with the same effective mass m, in agreement with the 
group velocity relation (38) which may also be writ¬ 


ten as dWjdv = vdP/dv. The corresponding energy- 
momentum relation reads 


W = W 0 + 


p2 

2m 


(41) 


at low momenta. Fig. [9] shows the numerically calculated 
energy-momentum relation and confirms the parabolic 
relation (411 at low momenta. In the opposite limit v 


v c the energy-momentum relation is linear 
W « v c P 


(42) 


which is also consistent with the group velocity relation 
(38 1 . Similar results have been recently obtained for a 


Q = 0 processing magnetic droplet in a model without 
DM interaction [33^ . 

In order to study the dynamics in the presence of an 
applied field we use Eq. (301 for the definition of the linear 
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FIG. 10: Linear momentum P as a function of time for a 
simulation where a field gradient (261 with g = — 0 . 001 , a = 
10 is applied to an initially static skyrmionium. The field is 
applied for 0 < t < 100 , and is then suddenly switched off. 




FIG. 11: Con tour plots for m 3 for a skyrmionium under the 
field-gradient (26 1 with g = —0.001, a = 10. (Left) The static 
Q = 0 skyrmionium of Fig. [3] is placed at the origin at initial 
time t = 0. (Right) The skyrmionium at simulation time t = 
160. It propagates along the x axis in the direction of the field 
gradient. The contour levels plotted are m 3 = 0.9, 0.6,0.3,0.0 
(solid lines) and m 3 = —0.3, —0.6, —0.9 (dashed lines). 


momentum and Eq. (231 applied for a field gradient (26 1 . 
We obtain 


P x = - J d x h( 1 - m 3 ) dxdy, 


Py = 0 


(43) 


which is the analogue of Newton’s law for the case of 
a non-topological (Q = 0) skyrmion. Under the field 
gradient the skyrmionium acquires a linear momentum 
along the gradient direction (x axis) according to Eq. (43) 


and confirms the Newtonian character of its dynamics. 
For a pure gradient h = gx we have 


Px = ~9H, 


Py = 0 


(44) 


where y = J (1 —m 3 ) dxdy is now the total magnetization 
of the skyrmionium. For the static skyrmionium of Fig. [3] 
we have y = 73 and thus P x = —73g. 

Finally, we have conducted a series of numerical simu¬ 
lations solving the initial value problem for the conserva¬ 
tive Landau-Lifshitz equation ©> using as initial condi¬ 
tion the static skyrmionium for n = 3. We apply a held 
gradient (26) with g = — 0 . 001 , a = 10 and find that 


the Q = 0 skyrmionium indeed propagates along the di¬ 
rection of the gradient (x axis) in sharp contrast to the 
Q = 1 skyrmion propagation along the y-axis shown in 
Fig. [ 6 ] 

The held gradient is applied for the time interval 
0 < t < 100 and is then switched oh. Fig. [To] shows 
the linear momentum as a function of time for the simu¬ 
lation. The increase is linear while the held is applied, in 
accordance with the prediction of Eq. (43) or (|44|. From 


the numerical data we obtain P x = 0.065 which remains 
almost constant throughout the simulation. This is in 
excellent agreement with the prediction obtained from 
Eq. (43) (deviation less that 1%). The prediction of the 
simplihed Eq. pTl) is R x = 0.073 and shows a deviation 


from the simulation result due to the approximation of h 
with a pure gradient. 

Fig. |TT] shows a contour plot for the magnetic con¬ 
figuration at the initial time t = 0 and at t = 160. We 
observe throughout the simulation that the configuration 
moves in a rather coherent way. It has evolved from the 
initial static skyrmionium (shown in Figs. [3] and [ 7 ] (left)) 
to a propagating one very similar to the those we calcu¬ 
lated earlier in this subsection (shown in Fig. 0 (right) 
for v = 0.07). 

At time t = 100 we suddenly switch off the applied 
field and observe a free motion of the skyrmionium. Its 
momentum remains constant and nonzero as shown in 
Fig. [TO] for t > 100. The skyrmionium propagates at a 
constant velocity confirming the Newtonian character of 
its dynamics. We switch off the held for various times 
and measure the energy, linear momentum and velocity 
of the freely moving skyrmionium. The results are shown 
by open circles which have been superimposed on the 
plots for the energy vs velocity and linear momentum vs 
velocity for steady-states in Fig. [9} They almost coincide 
with the data calculated for steady-state configurations. 


V. CONCLUSIONS 


We have given a theoretical description of the dynam¬ 
ics of topological and non-topological magnetic solitons 
in the presence of the Dzyaloshinskii-Moriya interaction. 
When a field-gradient is applied the Q = 1 skyrmion 
undergoes Hall motion perpendicular to the direction of 
the held gradient with a drift velocity proportional to the 
gradient, according to Eq. (27). When the applied held is 
switched off the guiding center ceases to move while the 
skyrmion undergoes cyclotron motion around its pinned 
guiding center. In contrast, the Q = 0 skyrmionium is 
accelerated in the direction of the gradient, with accel- 
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eration proportional to the applied field (see Eq. (431). 
When the applied field is switched off the skyrmionium 
continues traveling with constant velocity v in accordance 
with Newtonian dynamics of ordinary particles. 

The dramatically different dynamical behavior be¬ 
tween topological (Q ^ 0) and non-topological (Q = 0) 
solitons had been theoretically anticipated fTTriTTl but a 
full demonstration had not been given due to the lack 
of a system where both kinds of solitons could be conve¬ 
niently studied. The DM materials support robust topo¬ 
logical and non-topological magnetic solitons and thus 
appear to be particularly suitable for dynamical studies 
The dynamics can be induced by external fields, as in 
the present paper, or by magnon-skyrmion interactions 


The simulations presented show coherent propagation 
of solitons. However, significant distortions of the mag¬ 
netic configurations set in for larger held gradient. For 
the simulation of the Q = 0 skyrmionium we had to use 
a gradient as large as \g\ =0.1 in order to observe signif¬ 
icant distortions. 

We have also given a calculation of skyrmionium in 
steady-state propagation with velocities up to a maxi¬ 
mum critical velocity: v < v c . The energy-momentum 


relation resembles that of a non-relativistic particle for 
low momenta (see Eq. (411) but becomes relativistic-like 
(linear) for large morneta (see Eq. (421). 


The picture derived in this paper with theoretical tools 
together with straightforward numerical calculations is 
similar to the observed dynamics of magnetic bubbles [I] . 
Topological bubbles (Q ^ 0) are notorious for their skew 
deflection in the presence of an applied magnetic Held 
gradient. The deflection is in a direction almost perpen¬ 
dicular to the held gradient, in general agreement to the 
90° deflection of the Q = 1 skyrmion calculated in the 
present paper. The deviation from the 90° deflection is 
due to the presence of dissipation in bubble materials and 
should also be anticipated for DM skyrmions. Magnetic 
bubbles, stabilized by the magnetostatic interaction, are 
easily deformable under external probes while the long- 
range nature of the interaction makes numerical simu¬ 
lations more complicated. In contrast, chiral skyrmions 
are stabilized by the local DM interaction, thus they are 
robust and relatively easy to calculate numerically. The 
relation between chiral skyrmions and bubbles have been 
discussed in Refs. mmm, where the differences in the 
energetics, in stability and of their cores were stressed. 
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